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ABSTRACT 

Understanding the gas content of high redshift hales is crucial for studying the formation of the first 
generation of galaxies and reionization. Recently, Tseliakhovich & Hirata showed that the relative 
"stream" velocity between the dark matter and baryons at the time of recombination - formally a 
second order effect, but an unusually large one - can influence the later structure formation history of 
the Universe. We quantify the effect of the stream velocity on the so-called "characteristic mass" - the 
minimum mass of a dark matter halo capable of retaining most of its baryons throughout its formation 
epoch - using three different high-resolution sets of cosmological simulations (with separate transfer 
functions for baryons and dark matter) that vary in box size, particle number, and the value of the 
relative velocity between the dark matter and baryons. In order to understand this effect theoretically, 
we generalize the linear theory filtering mass to properly account for the difference between the dark 
matter and baryonic density fluctuation evolution induced by the stream velocity. We show that the 
new filtering mass provides an accurate estimate for the characteristic mass, while other theoretical 
ansatzes for the characteristic mass are substantially less precise. 



1. INTRODUCTION 

Gas rich dark matter halos in the early universe 
serve as a nu rturing ground for dwarf galaxi es (e.g. , 
Ricotti et all l2002bl lal: iBromm et aLl l2002l [l999l: 



Abel et al 



.2002t iNaoz et al.l 120061 lYoshida et al I 
20061 120081: IGreif et al.l 120101: IClark et all 120111 
Bromm fc Yoshidal 1201 IL and references therein). 



Their properties are important to quantify, as they are 



at the onset of structure formation fe.E.. iShaoiro et al. 


2004 


; iCiardi et al.ll2006 


; iHoeft et al.ll2006l: Gnedin et al. 


2008 


: lOkamoto et al.l 


20081: iTrenti & Stiavellil 120091). 



More than that, even if the smallest of gas rich halos 
are too small for efficiently cooling via atomic hy- 
drogen lines and may not host actual galaxies, these 
"mini-halos" may p roduce a 21-cm sign a ture in future 
radio o bservations (iKuhlen et al.l (I2006D: iShapiro et al. 
'2006h : iNao^ fc Barkanal (j2008li " but see lFurlanetto fc OhI 



2006!)) and might block some of the ionizing radiation, 
causing an ov erall delay in the initia l progress of reion- 
izatio n fe.g.. iBarkana fc Loebl 120021: llliev et al.l l2003al 
120051 iMcOuinn et al.l I2007D . Thus, the evolution of 
the gas fraction of dark matter halos at various epochs 
during the early evolution of the universe is of prime 

importance.^ 

Recently, iTseliakhovich fc Hiratal ()2010f ) showed that 
not only the amplitudes of the dark matter and bary- 
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onic density fluctuations were different at early times, 
but also were their velocities. After recombination, the 
sound speed of the baryons dropped dramatically, while 
the dark matter velocity remained high - thus, the rela- 
tive velocity of bary ons with respect to the dark m atter 
became supersonic. ITseliakhovich fc Hiratal ()2010D also 
showed that this relative velocity between the baryons 
and the dark matter remained coherent on scales of a 
few mega-parsec and was of the order of ^ 30 km scc~^ 
at the time of recombination. This relative velocity is 
often called the "stream velocity" in the literature, and 
throughout this paper we will use this term. The stream 
velocity effect has previously been overlooked, because 
the velocity terms are formally of the second order in 
the perturbation theory and should be neglected in the 
linear approximation. However, this second order ef- 
fect is unusually large, resulting in the numerically non- 
negligible suppression of power at mass scales that cor- 
respond to the first b ound objects in the Universe (e.g., 

lYoshida et a"Lll2003al) . 

Using the Prcss-Schcchtcr (Press fc S chechteri [TqtI 
formalism, [Tse liakhovic h fc Hirata (20101 ) showed that 
the number density of halos is reduced by more than 
60% for halos wit h M = 10^ Mf^ at z = 40. In a 
subsequent paper, ITseliakhovich et all ()2010[ ) also in- 
cluded the baryonic tem perature fluctuations following 
iNaoz fc Barkanal (|2005[ ). They found that the stream 
velocity also resulted in much higher "characteristic" 
mass - the minimum mass for a dark matter halo ca- 
pable of retaining most of its gas - as compared to the 
case without the stream velocity (e.g. iNaoz fc Barkanal 
|2007|) . As has been shown in subsequent studies, the 
stream velocity eff ect has imp ortant implicatio ns on 
the first structur e s (IStacv et al.lf 2011: Maio et al. 
Greif et al.l I2OIII: INaoz et all 1201 2: Fialkov c tlT 



O'Learv fc McOuinnl l2012t iBovv fc Dvorkin 120121 ) and 



2011 



2011 



may also affect the redshifted cosmolog i cal 21-cm signal 



(iDalal et al.ll2010t iBittncr fc Loebll201lL 
IVisbal et a'Lll20ia IMcOuinn fc O'Learv! 



Yoo et al.ll20nl: 
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In this paper we explore the effect of stream velocity 
on the gas fraction in dark matter halos and compare 
the simulation results to the predic tions from the lin- 
ear the ory (e.g. iTseliakhovich et alJ 120101. In our first 
paper (|Naoz et alj 120121 hereafter Paper I) we quanti- 
fied the stream velocity effect on the evolution of the 
halo mass function with cosmological simulations. We 
used three different sets of high resolution simulations in 
order to study the stream velocity effect systematically, 
thus understanding the overall trends (instead of concen- 
trating on specific halos). We used a set of simulations 
with different box sizes, particle numbers, and the values 
for the stream velocity to analyze the suppression of the 
structure formation as a function of the stream veloc- 
ity. In Paper I we found that the total number density 
of halos is suppressed by ^ 20% at z = 25 in regions 
of the universe that happen to have Vbc = Icvbc, where 
Uvbc is the (scale independent) rms fluctuation of the 
stream velocity on small scales. In rare patches where 
Vhc = 3.4(Tvbc, the relative suppression at the same red- 
shift reaches 50%, remaining at or above the 30% level 
all the way to z = 11. Perhaps the most interesting 
phenomenon that we found was the high abundance of 
"empty halos", i.e., halos that had their gas fra_ctions be- 
low half of the cosmic mean baryonic fraction /b. Specif- 
ically, we found that for Ubc ~ Icvbc cill halos below 
10^ Mq are empty aX z > 19. As a result, the high abun- 
dance of empty halos can significantly delay the forma- 
tion of gas rich "minihalos" and the first galaxies. In this 
paper we investigate the effect of the stream velocity on 
the gas fraction in halos. In particular, we quantify the 
dependence of the characteristic mass on the magnitude 
of the stream velocity. 

For completeness we first describe the parameters and 
initial conditions of our simulations in fj2] We present 
our results and analysis of the gas fraction in halos and 
comparison to the linear approximation in Section |31 Fi- 
nally we offer a brief discussion in gl 

Throughout this paper, we adopt the following cosmo- 
logical parameters: (f^A, ^^M, ^b, n, ggi Hq)= (0.72, 0.28 , 
0.046, 1, 0.82, 70 km s"! Mpc"!) ()Komatsu et al.|[2009l ). 

2. THE SIMULATIONS 

2.1. Basic Parameters and Settings 

In this work we use a parallel iV-body/hydrodynamic s 
code GADGET-2 (|Springel et all 120011 ISpringell l2005h . 
Below we describe the general features of our 3 simulation 
sets, which are also summarized in table [T] 

1. The first set, named 'W = 256", uses a total of 
2 X 256^ dark matter and gas particles within a 
cubic box of 200 comoving kpc on a side. To real- 
ize statistically significant number of halos in such 
a small box, we artificially increase gravitational 
clustering in the simulation by setting erg = 1.4. 
We choose this box size so that a lO'' halo is 
resolved with ~ 500 particles - the valu e needed to 
estimate the halo gas fraction reliably (jNaoz et al.l 
I2009D . The gravitational softening is set to be 
40 comoving pc, well below the virial radius of a 
10^ M0 halo (~ 680 comoving pc). All the simula- 
tions in this set are initialized ai z = 199. 

2. The second set, named 'W = 512", uses a total 
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Parameters of the simulations 
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(Tvbc stream 




1 1 

km sec 


velocity 


256 runs, 


0.2 Mpc, 


^in — 199 


256o 








256io-+p 


5.8 


1 


256i.7<T 


10 


1.7 


2563.4,7 


20 


3.4 


512 runs, 


0.7 Mpc, 


^in — 199 


512o 








5121,, 


5.8 


1 


512i.7,T 


10 


1.7 


5123.4,T 


20 


3.4 


768 runs, 


2 Mpc, 


•^in — 99 


768o 








768ia 


3 


1 


768 1.7a 


5 


1.7 


7683.4a 


10 


3.4 



of 2 X 512^ dark matter and gas particles within 
a cubic box with the size of 700 kpc. In this set 
we also artificially increase erg to 1.4. With these 
parameters, a halo with 500 dark matter particles 
has a mass of 5 x 10'^ M0. The softening length 
is set to be 68 comoving pc. All the simulations in 
this set are initialized at z = 199. 

3. The final set of simulations uses 2 x 768"^ dark mat- 
ter and gas particles (which we name the 'W ~ 
768" set) in a 2 Mpc box, and starts at z = 99. 
For these parameters a halo with 500 dark matter 
particle has a mass of ~ 10^ Mq. The softening 
length is set to be 0.2 comoving kpc. We use the 
"correct" value of erg = 0.82 for this simulation set. 

In each simulation set, we explore a range of the values 
for the stream velocity (see table [1]). 

2.2. Initial Conditions 

As has been shown by iNaoz etldl ()2006f ) and 
iNaoz k, Barkanal (|2007[) , setting up initial conditions for 
cosmological simulations on small spatial scales is a del- 
icate issue. High accuracy in initial conditions is crucial 
for accurately predicting the halo mass function in the 
lowest mass regime (M < 10^ M fri). 

Following INaoz et al.l (|2011[ ). we generate separate 
transfer fu nctions for dark matter and baryons as de- 
scribed in iNaoz fc Barkanal ()2005[ ). Ideally, streaming 
velocities should be realized in the initial conditions in 
a self-consistent way with the transfer functions that are 
calculated up to the second order in the perturbation 
theory (c.f. recent studies bv lO'Learv fc McQuinnll2012l : 
iMcQuinn fc O'Learvl2012|) . However, that would require 
computing transfer functions up to the second order in 
perturbation theory, and those are not readily available. 
Instead, similar to all previous simulation studies, we 
used the transfer functions computed in the linear ap- 
proximation. We account for that choice in i j3.3| when 
we compare our simulations to the perturbations theory, 
and thus our results remain self-consistent. 

For all runs, glass-like initial conditions were gener- 
ated using Zel'dovich approximation. For baryons, we 
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have used a glass file with positions shifted by a random 
vector, thus removing artifici al coupling between n earby 
dark matter and gas particles (lYoshida et al.|[20"03bf ). We 
note that we have used the same phases for dark matter 
and baryons in all of our simulations, since we showed in 
Paper I that the spatial shift between baryons and dark 
matter is unimportant. 

iTseliakhovich fc Hiratal (|201Clf ) demonstrated that, 
while the stream velocity varies in space, its coherence 
length is quite large, many Mpc. Hence, on scale of our 
simulation boxes, it can be treated as constant bulk mo- 
tion of baryons with respect to the dark matter. We 
include the effect of stream velocity by adding, at the 
initial redshift, an additional velocity to the x compo- 
nent of the baryons velocity vector. We test a range 
of values for the stream velocity, which is convenient to 
quantify in terms of its rms value on small scales, a^bc- 
Specifically, we test Vhc = Icvbc through Wbc = 3.4CTvbc 
for all the simulations sets (see table [1]). 

2.3. Halo Definition 

We locate dark matter halos by running a friends-of- 
friends group finder algorithm with a linking parameter 
of 0.2 (only for the dark matter component). We use 
the identified particle groups to find the center of mass 
of each halo. After the center is located, we calculate 
density profiles of dark matter and baryons separately, 
assuming a spherical halo and using 2000 radial bins be- 
tween r,iii,i = kpc and Tmax = 20 kpc. Using the den- 
sity profiles, we find the virial radius r^ir at which the 
total overdensity is 200 times the mean background den- 
sity, and compute the mass and th e gas fraction of each 
halo w ithin that radius. Recently, lO'Learv fc McQuinnI 
()2012|) used an unconventional definition for halos by us- 
ing the highest baryonic density peaks as the center their 
halos. This method preferentially results in larger gas 
fraction for (Jvbc 7^ 0, compare to our more conservative 
method. 

Recently, IMore et all (j2011[ ) showed that halos identi- 
fied by the friends-of-friends algorithm enclose an average 
overdensity that is substantially larger than 200, and its 
specific value depends on the halo concentration. In our 
approach we use the friends-of-friends algorithm only to 
find the center of mass of a halo, and compute the actual 
halo mass using the spherical overdensity of 200. 

We only retain halos that contain at least 500 dark 
matter particles within their virial radii. The choice al- 
lows us to estimate halo masses to about 15% precision 
(jTrenti et al.l I2010D and to estimate halo gas fraction s 
reliably to a similar level of accuracy (jNaoz eTall [20091) . 
However, for some of our fit calculations we also include 
halos with the number of particles as low as 100; if we 
do that, we assign a lower weight in the fit t o these ha- 
los ac cording to the resolution study done in iNaoz et"al] 
(|2009l ). see Appendix [X] for more details. 

3. RESULTS 

Let us consider the various scales involved in the for- 
mation of cosmic structure. On large scales gravity dom- 
inates other forces and gas pressure can be neglected. 
On small scales, on the other hand, the pressure dom- 
inates gravity and prevents baryon density fluctuations 
from growing together with the dark matter fluctuations. 



The relative force balance at a giv en time can be char- 
acterized by the "filtering scale" (|Gnedin fc Hull [19981 ) 
- a physical scale above which a small gas perturbation 
can grow due to gravity overcoming the pressure gradi- 
ent. In the non-evolving background, the filtering scale 
coincides with the classical [jcanSi (,1928[) scale, but in the 
expanding universe the two scales typically differ by a 
significant factor. 

Immediately after recombination Compton scattering 
of Cosmic Microwave Background (CMB) photons on the 
residual free electrons after cosmic recombination kept 
the gas temperature coupled to that of the CMB, the 
Jeans mass was constant in time and equal to the fil- 
tering scale. However, at z 130, the gas temperature 
decoupled from the CMB temperature, the Jeans mass 
began to decrease with time as the gas cooled adiabat- 
ically, and the filtering scale lagged behind the Jeans 
scale. 

Based on re sults from an early numerical simulations, 
iGnedinI (|2000[ ) suggested that the filtering mass also de- 
scribes the largest halo mass whose gas content is sig- 
nificantly suppressed compared to the cosmic baryon 
fraction. The latter mass scale, commonly called the 
"characteristic mass", is defined as the halo mass for 
which the enclosed baryon fraction equals half of the 
cosmic mean. Thus, the characteristic mass distin- 
guishes between gas-rich and gas-poor halos. Many semi- 
analytical models of dwarfs galaxies use the characteris- 
tic mass sc ale in order to estimate the gas fraction in 
halos (e.g IBullock et al.l [2000l: IBenson et al.l [2002a|[bl: 
iSomervillel 120021: iBovv fc Dvorkinll2012[) . Theoreticallv 
this sets an approximate minimum value on the mass 
that can still form stars. 

3.1. Non-linear Behavior: the Characteristic Mass 

For halos, IGnedinI ()2000f ) defined a characteristic mass 
Mc for which a halo contains half the mean cosmic baryon 
fraction /f,. In his simulation he found the mean gas 
fraction in halos of a given total mass M, and fitted the 
simulation results to the following formula: 



fg,ca\c fh.O 



1 + (2"/3 _ 1^ 



IT 



-3/q 



(1) 



where /b,o is the gas fraction in the high-mass limit. In 
this function, a higher a causes a sharper transition be- 
tween the high-mass (constant fg) limit and the low-mass 
hmit (assumed to be fg a M^). IGnedinI (120001 ) found a 
good fit for a = 1, with a characteristic mass that in fact 
equaled the filtering mass by his definition. 

The characteristic mass is essentially a non-linear ver- 
sion of the filtering mass, and so it also measures the com- 
petition between gravity and pressure. At high masses, 
where pressure is unimportant, fg — )■ /b,0: while the low 
mass tail is determined by the suppressio n of gas ac- 
cretion by gas pressure. INaoz et alT (|2009[ ) found that 
the filtering mass from linear theory (calculated in a self 
consistent way) is consistent with the characteristic mass 
fitted from the simulations, for two (pre-reionization) sce- 
narios that they tested: with no stellar heating and 
a case of a sudden fiash of stellar heating at a given red- 
shift. In a foUowup paper. [Naoz et ahl ( 20111) found the 
same agreement between the the linear and non-linear 
theory, and showed that alternative initial conditions 
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Fig. 1. — Binned gas fraction in the "N=512" runs at redshifts 25, 
19, and 15 (from bottom to top panels) for 4 values of the stream 
velocity j^bc = 0, v^c = lo-vbc, fbc = l-7crvbc, and iJbc = 3.4(Tvbc 
(blue squares, purple diamonds, brown stars and red triangles, re- 
spectively). We also show the best fit as evaluated from equations 
fill and |2J, with dotted and solid lines respectively. (See also Fig- 
ure[7]for the gas fraction behavior for all halos, i.e., not binned.) 

models yield a different (higher by about 50 per cent) 
minimum mass (both the linear and non- linear), since 
the system retains a memory of the initial conditions. 
For clarity, we emphasize that the statement {Mc = Mp) 
refers to our definition of Mp in equation (jS]). 

In Figure [T] we present the gas fraction as a function of 
halo mass for our "N=512" simulation selQ. As one can 
see, the halo gas fraction drops dramatically at lower halo 
masses for large values of the stream velocity. This trend 
introduces a qualitatively different behavior in the gas 
fraction as a function of halo mass, which is not captured 
by equation ([1]). The best fit ansatz [Eq. ([T])] is shown 
with dotted lines in Fig.[TJ It clearly does not capture the 
behavior of the gas fraction as a function of halo mass for 
'^bc 7^ 0. Therefore we introduce a new fitting formula 
for the gas fraction as a function of halo mass, 



fg,ca.lc fh.O 



1 + (2^ - 1) 



M, 
IT 



-1/7 



(2) 



The new fitting formula reduces to Equation ([T]) for 7 = 
3/3 = a. 

Although we add another free parameter to the fit, 
finding the best fit model presents some hurdles which we 
discuss in the Appendix. The best-fit value of from 
Equation ([2]) is the same as the value from Equation ([T]) 
for the Vbc = case, as can be expected from the fact that 
Equation ([T]) provides a good fit to the simulation results 
in the Vhc = case. It is interesting, however, that the 
best fit values for Mc from Equation ([2]) and Equation 
([T|) are still the same even for large values of Wbc, as can 

In Paper I we showed that, qualitatively, all of our simulation 
sets behave similarly as a function of mass and redshift. There- 
fore, to avoid redundancy, we show here the gas fractions for only 
"N=512" runs. 



be seen in Figure [S] However, the new fitting formula 
gives a better overall fit for high Ubc values, especially in 
the high mass limit. 
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Fig. 2. — Best fit values for the characteristic mass and limiting 
baryon fraction /t, g as functions of redshift for various values of 
the stream velocity d^c = Oi ''be = l^vbci ''be = l-7o"vbci and 
"be = 3.4(Tvbe (blue squares, purple diamonds, brown stars and 
red triangles respectively) for our A' = 762 simulation set. We 
also show the evolution of the fully self-consistent filtering mass 
(Equation [6} with solid curves. 




Fig. 3. — Same as Figure [2l but for the A' = 512 simulation set. 

Figures [JHH show the best- fit values for the character- 
istic mass and limiting baryon fraction /b,o at a range 
of redshifts. At the highest redshifts and large values of 
(Tvbc rnost of the halos are empty halos (i.e., halos with 
gas fraction lower then the half of the mean cosmic bary- 
onic fraction, see Paper I, figures 2-4). Therefore, no 
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perature {St)' 



Fig. 4. — Same as Figure |2] but for the N = 256 simulation set. 



convergence could be achieved in the fitting procedure 
and the parameters of the fit cannot be measured. Fur- 
thermore, as has been noted in Paper I, the N = 256 set 
suffers from poor statistics for M > IO^Mq, resulting in 
large error bars for the best-fit values of the parameters. 

An important point to make is that at low redshifts the 
baryon fraction at the highest mass bins (for all runs) ap- 
proaches the same value irrespectively of the magnitude 
of the stream velocity. This is, of course, expected, as 
the global trend for the baryon fraction is to approach 
the cosmic mean (see Appendix 1X1 Figure [9]). In fact, 
in the case of the largest stream velocity we consider, 
Vhc = 3.4(Tvbc, the N = 256 and N = 512 simulation 
sets do not quite reach the cosmic mean values expected. 
This is most likely due to the low abundance of the most 
massive halos, as can be seen in Figure 4 of Paper I. The 
N = 768 simulation set, however, does not suffer from 
that incomplete convergence. 

3.2. Linear Theory Predictions: the Filtering Mass 

In the li near approximation, the filtering mass, first 
defined bv IGnedin fc HuH (jl998D . describes the highest 
mass scale on which the baryon density fiuctuations are 
suppressed si gnificantly compared to the dark matter 
fiuctuations. IGnedin fc Hull ()1998D only considered the 
low redshift case, where the baryonic and dark matter 
fiuctuations have the same amplitude at large scales. 
INaoz fc Barkanal (|2007D relaxed that assumption and ex- 
tended the computation of the filtering mass to early 
times, during which the amplitude of the baryonic fiuc- 
tuations is below the amplitude of the dark matter fiuctu- 
ations even on large scales. Both studies, however, only 
considered a case of zero stream velocity. 

In order to extend the derivation of the filtering mass 
to the non-zero stream velocity case, we first introduce 
the coupled second order differential equations that gov- 
ern the evolution of the density fiuctuations of the dark 
matter (^dm), and the baryons ((5b) and the baryon tem- 



2i 



Sdra + 2i/(5dm — /dm — Vbc ' k^dm = 

a 

3 VL. 

2^0 (./b<5b + /dm^dm) + 



Vbc • k 



^dn 



Sb + 2H5b 



(3) 



(4) 



-i/n — !^ {fhSh + /dm^dm) — ' , 

2 a'' 



where Vim is the present day matter density as a fraction 
of the critical density, k is the comoving wavcnumber, a 
is the scale factor, /i is the mean molecular weight, Hq is 
the present day value of the Hubble parameter H , and T 
and 5t are the mean baryon temperature and its dimen- 
sionless fiuctuation, respectively . These equations are a 
compa ct form of equations 5 in iTseliakhovich fc Hiratal 
(|20Ta) . where we used the fact that Wbc oc 1/a, and the 
baryon equation includes the pressure term whose form 
comes from the equation of state of an ideal gas. The 
linear evolution of the tem perature fiuctuati ons is given 
by (jBarkana fc Loebl [20051: iNaoz fc Barkanall2005D 



(ISt 



2 dS^ 



Xe{t) _4 
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where Xe{t) is the free electron fraction as a function 



-y is the photon density fluctuation, 
^ and and St^ arc the mean 



of cosmic time t, S. 
t^ = 8.55 X 10-i3yr 

photon temperature and its dimensionless fluctuation, 
respectively. Equation ([5|) describes the evolution of the 
gas temperature in the post-recombination era, but be- 
fore formation of first galaxies, when the only external 
heating arises from Compton scattering of the remaining 
free electrons on the CMB photons. The first term in 
Equation ([5]) comes from the adiabatic cooling or heat- 
ing of the gas, while the second term is the result of the 
Compton interaction. 

In the top right panel of Figure [5] we show an example 
of the solution of Equations We plot the ratio for 
Sh/Stot as a function of the wavenumber k for the fully 
self-consistent linear calculation, i.e., starting at the time 
of re combination and using the exa ct transfer functions 
from ITseliakhovich fc Hiratal ()2010[ ) . We consider cases 
with Wbc = 0, Vhc = Cvbc, and Vhc = 2(Tvbc at z = 10 and 
z = 25. For larger Vhc values, the drop in Sh/Stot occurs 
at larger scales, i.e., the suppression of the baryonic per- 
turbations relative to the total matter fiuctuations shifts 
to larger m asses. 

Following INaoz fc Barkanal ()20Q7D . we re-define the fil- 
tering scale (specifically, the filtering wavenumber kp) to 
include the stream velocity effect as 



Stot 



1 + rLss 



(6) 



where v = 7;bc/cvbc and (ivbc is the (scale-independent) 
rms of the stream velocity at small scales. The parameter 
fLSS (a negative quantity) describes t he relative differ- 
ence between Sh and Stot on large scales (|Naoz fc Barkan"al 
[2003), i.e., 

?'LSS = T , (7) 

"tot 
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Fig. 5. — Fully self-con sistent linear calculation of the growth of the dark matter and baryonic fluctuations (i.e., following 
ITseliakhovich fc Hiratall2010l ). We show in the bottom panel the filtering mass as a function of redshift in the regions with = 0, 1 and 
2 (Tvbc (see labels). In the top panels we show 5b/^tot as a function of the wave number k. In the right top panel we plot (5b/'5tot in the 
regions with v^,^ = (dotted red lines), IfTvbc (solid magenta lines) and 2(t^,i,c (dashed blue lines) at 2: = 10 (top lines) and z = 25 (bottom 
lines), see labels. In the top left panel we show an example for the d^c = IcTy bc c ase at 2 = 15 and plot the two fitting functions (see text 
for details), one which reproduces the drop of i5b/<5tot with wavenumber, Eq. II14I I (red dashed line) and the fitting function to the second 
order, i.e., Eq. l(6]l, (brown short dashed line). 



where A = (5b — 6tot (see also iBarkana fc Loebl 120051 ) . 
The fihcriiig mass is defined from kp simply as: 



47r _ /I 27r 



(8) 



where po is the mean matter density today. 
To find kp in a general case, we write it in the form 



fc2 



1 



pit) = ^ 

u{t) 



(9) 



where u(t) is to be determined. Then, using equation 
([6]), we expand the baryonic fluctuation as a function of 
wavenumber k. 



A 



LSS 



0{k^), 



(10) 



where Alss = ^Lss^tot [eq- (0)] obeys the following equa- 
tion to the first order of fc, 



A 



LSS 



2HA 



LSS 



— /dn 

a 



iVbc • k(5d„ 



(11) 



Note that in the case of Wbc = 0, the linear term of fc 
has a zero coefficient, and thus the right hand side of 



this equation is simply zero (see IBarkana fc Loebl 120051: 
iNaoz fc Barkanall2007n . Substituting the expansion from 
equation ([TU| into equation ([5]), and using equations ^ 
and (jlip . we obtain an equation for u: 



2i/u=/dn.(l 



1 fceT 



((5b + 6t) (12) 



Vbc 



In the limit of Ubc = (i.e., ly — 0) thi s equat ion reduces 
to equation (12) of INaoz fc Barkaiial (|2007f ). and thus 
results in the same filtering mass found in that study. 
We can solve Equation ([TS]) to find u{t), 



u{t) = 

/dm (1 



dt" 



dt 



,kBf{t') 



(13) 

6T{t')) 



dt" 



l(t") 



dt' 



tro 



Vbc 



(i')-k) (5dm(i') 



where Vbc,rec is the stream velocity at the moment of 
recombination t^ec, so that Vbc(i) = Vbc,reca(^rec)/a(i)i 
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and k — k/fc is the unit wavenumber vector. In the 
bottom panel of Figure [5] we show the evolution of the 
filtering mass as a function of redshift for t;bc = 0, 1 and 
2(Tvbc- The values of the filtering mass for Ubc 7^ as de- 
fined by Equation ^ are larger up to an order of mag- 
ni tude at high redshif t s as c ompared to the definition 
of iTseliakhovich et al.l ()2010[) . We emphasize that this 
difference is entirely due to the diffe rent definition of the 
filterin g scale, not due to any error in lTseliakhovich et alj 
(|2010( ) calculations. 

The filtering scale kp can also be obtained simply by 
fitting equation ([6]) to the calculated values of Jdm and Sb, 
using equations ©-(E]). iNaoz fc Barkanal (I2007D found 
a functional form that can be used to produce a good fit 
for the drop of the wavenumber. Generalizing it to the 
case of stream velocity we write: 



f .(1 

Otot 



rhss) 



1 

1 + - 



1 



nl + TLss 1 + 



(14) 



and n must be adjusted at each redshift. In the top left 
panel of Figure [5] we compare this fitting formula to the 
fully self-consistent linear calculation (long dashed line) 
for which reproduce the drop of (5b/(Stot as a function of 
k fairly well. For the example considered in the figure, 
i.e., z ~ 15 and v^c = Ifvbc; we find n = 0.46 and 
kp = 253.9 Mpc""'^. We also show the resulted fit using 
the second order in k approximation, i.e. Eq. ([6|). 

3.3. Comparison Between the Linear Theory 
Predictions and the Nonlinear Results 

In order to compare the filtering mass to the char- 
acteristic mass we calculated the fi ltering mas s in a 
self c onsistent way, as was done in iNaoz et all (|2009L 
2011f). In other words we use the transfer function from 



Naoz fc Barkanal (|2005[ ) with a boosted velocity for the 
baryons at z = 99 (z = 199) for the N = 768 {N = 512 
and = 256) set as initial conditions. We then evolve 
the dark matter and baryon in time according to equa- 
tions ©-(IS]). We note that in all our calculations we 
included the fact that the boost of the velocity was in- 
cluded in the simulation only in one axis, thus terms 
which are proportional to Vbc • k are reduced by a factor 
3 compare to the global average. 

Our simulation sets N = 256 and N ~ 512 are 
initialized at redshift 199 at which Compton heating 
by the CMB photons significantly affects the evolution 
of the linear modes and spec ifically the filtering mass 
()Naoz fc Barkanal [20051 [2001 . However, GADGET-2 
does not include CMB Compton heating. Hence, to 
compare apples and apples, we neglected the Compton 
heating contribution to the filtering mass when compar- 
ing these two simulation sets to the linear approxima- 
tion. This is the reason that our values of the filtering 
mass in Figures [3| and [Jj are lower than the values for the 
no str eam velocity case in Figure 3 of INaoz &: Barkanal 
([2OOI . 

We show the linear theory filtering mass as a function 
of redshift for all the cases we consider in Figures [2H11 
The filtering mass is consistent with the characteristic 
mass within our fit errors for all simu lation sets. There- 
fore, we conclude in agreement with INaoz et ahl (120091 
l2011t) , that the evolution of the characteristic mass can 
be understood using the linear approximation predictions 
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Fig. 6. — Best fit values for the characteristic mass and limiting 
baryon fraction /b,o as functions of redshift (same as in Fig. [Ill, but 
now compared with previou s theoretical mode l s for the character- 
istic mass. Solid lines sh ow INaoz fc Barkanal 1120071 ') definition of 
the filtering mass (Eq. I15|l that neglects the stream velocity contri- 
bution. We also plot M^j j (dashed lines) and Mesc (dotted lines) 
as alternative models, see text for details. 

for the filtering mass. 

3.4. Comparison with the Previous Definition of the 
Filtering Mass and Other Mass Scales 

Recently ITseliakhovich et al.l (|2010( ) showed that in- 
cluding the effects of the relative velocity between the 
dark matter and the baryons at recombinations results 
in a higher filtering mass as compared to the case of 
'^bc = (by about o rder of magnitude for the global av- 
erage). INaoz fc Bark ana (2007) defined the filtering scale 
in the case of Ubc = as: 



1- — 

Stot k% 



■ fLSS 



(15) 



In Figure [6] we show thus defined filtering mass for 
the N = 512 s i mulat ion set. As can be seen, the 
INaoz fc Barkanal (|2007f ) definition underestimates the 
characteristic mass in the high Vhc limit, since it neglects 
the difference between the dark matter and baryons den- 
sity fluctua tions. 

Recentlv lStacv et al.l (|201lD and lNaiman et al.l (|201lD 
suggested that, given a high initial stream velocity, the 
baryon evolution is dominated by the relative motion 
of dark matter and gas, and thus the gas sound speed 
should be replaced with the effective sound speed Wcff , 



(16) 



where Wbc(2) = Wbc, 0/(1 + 2), in the Jeans mass definition. 
This effective Jeans scale kj^s can be written as 



(17) 



where G is the gravitational constant and pm is the av- 
erage density. The effective Jean mass A/eff associated 
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with this scale length is simply 



47r _ / 1 27r 
Mcff = —po T^T 



(18) 



We show this mass scale in figures IH] with dashed lines. 
As one can see, this mass scale overestimates the char- 
acteristic mass at all redshifts and for all values of the 
stream velocity that we considered. This is not unex- 
pected, since the Jeans mass always overestimates the 
scale at which pressure starts to overtake gr avity in the 
expanding background (jGnedin fc Hui|[l998l ). 

The evolution of the characteristic mass as a function 
of redshift can also be modeled by considering the escape 
velocity of the gas. Given a mass of a halo M, the escape 
velocity is simply Vesc = \J^GMjr^ where r is the virial 
radius of the halo. For example, for a halo of 10^ the 
escape velocity is about 0.77 km sec~^, while the stream 
velocity for Wbc = 3.4crvbc at z = 15 is 1.6 km sec~^ and 
at z = 25 is 2.6 km sec~^. Thus, it is not surprising that 
halos below 10^ are empty in that redshift range - 
the stream velocity is simply much larger then the halo 
escape velocity, so the dark matter halo is unable to ac- 
crete any gas. We can, thus, estimate a halo mass Mosc 
below which the stream velocity is larger then the escape 
velocity. 



v/(2GHo)2a.A, ' 



(19) 



where Ac = 200 is the virial overdensity and v\,c(z) is 
the stream velocity at redshift z. We show this limit 
in Figure IH] with dotted lines. Equation ((T^ predicts a 
much stronger evolution of the characteristic mass than 
is actually observed in our simulations and provides a 
poor fit to simulations results. 

4. CONCLUSIONS 

We have used three-dimensional hydrodynamical sim- 
ulations to investigate the effects of stream velocity on 
the gas f raction in h i gh re dshift halos. In a compan- 
ion paper iNaoz et all ([20121) . we studied the effect of the 
stream velocity on the total halo mass function. In this 
work we focus on the effect of the stream velocity on the 
gas fraction in halos and on the evolution of the charac- 
teristic mass, and compare the simulation results to the 
linear approximation. 

In a first improvement over the earlier results, we intro- 
duce a new fitting formula (Eq. [2]) which offers a much 
better fit to the gas fraction as a function of halo mass 
at a given redshift in the limit of large stream velocities, 
while returning essentially the same values of the charac- 
teristic mass Mc as the previously used functional form 
(see Figures [U [7] an d[8l and see Appendix IX)) . 

Previous studies (jNaoz et al.ll2009L l201lf rshowed that 
a quantit y defined in the linea r approximation, the filter- 
ing scale (|Gnedin fc Huil[T998h . provides a good match to 
the nonlinear characteristic mass measured in numerical 
simulations. We introduce a new definition for the lin- 
ear filtering mass that accounts for two effects neglected 



in lGnedin fc Hull (I1998D : the deviation of the amplitude 
of baryonic fluctuations from the dar k matter fluctua- 
tions on large scales (considered flrst bv lNaoz fc Barkanal 
l2007f l and the stream velocity between the dark matter 
and baryons on small scales, which we include in the defi- 
nition of the filtering mass for the first time in this paper. 
The latter effect may result in the filtering mass being 
up to an order of magnitude larger at high redshifts for 
high values of the stream velocity, as compared to the 
case when the stream velocity is neglected. 

Finally, in comparing our simulations results to the lin- 
ear calculation (using our new definition of the filtering 
mass) , we find that the filtering mass (i.e. a linear quan- 
tity) offers an accurate match to the actual nonlinear 
characteristic mass measured from the simulations, at all 
redshifts and for all values of the stream velocity that we 
simulated. On the contrary, previous theoretical models 
that used as the characteristic mass scale either the halo 
mass with the escape velocity equal to the stream veloc- 
ity or the Jeans mass for the "effective" gas sound speed 
provide only poor fits to the simulation results. 

It has been suggested in the literature that gas rich 
low mass halos may play an important role in cos- 
mic reionization, and that they can p r oduce distinct 
21-crn signatures (iKuhlen et al l ()2006l): IShapiro et al. 
' 20061) ; INaoz fc Barkanal (|2008[) but see lFurlanetto fc Oh 
20061) ). For example, minihalos (halos of mass ^ 



10 Mq) can potentially block ionizing radiation and in- 
duce an over all delay in the initial p rogress of reion- 
ization (e.g.. IShapiro fc Girouxl 119871: Barkana fc Loeb 

20041 : llliev eTal 



2002 




2005 





Iliev et al. 200 3b: Sh apiro et al.l 



McOuinn et al.ll2007l ). However, our results here 



suggest that at high redshifts the stream velocity effect 
results in large variations in the characteristic mass - i.e. 
the minimum mass of a gas ric h halo. Thus, if re ion- 
ization started sufficiently early (jYoshida et al.|[2"007l) , in 
patches of the universe where the stream velocity is large 
there were fewer gas rich halos that can absorb ionizing 
photons. Hence, in these patches the delay of the reion- 
ization caused by minihalos would be less than in regions 
that happen to have a small value of the stream velocity 
and, hence, a large abundance of minihalos. Therefore, 
not only the formation of the first generations of galax- 
ies may be affected by the stream velocity effect, but also 
the whole process of reionization may proceed differently 
in regions with very different strea m velocities. This ef - 
fect has been considered recentl y bv lVisbal et al.l (|2012[ ) 
and lMcQuinn fc O'Learvl (|2012D . but our results indicate 
that it can even be stronger than previously estimated. 
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APPENDIX 
FIT CALCULATION 

As can see in Figure [U eq. ([1] dotted lines in that Figure) dose not capture the behavior of the gas fraction as a 
function of mass. During our search of better fit and a better formula for the fit we found that the low mass tail has a 
significant effect on the fit for large Vhc values. It is not surprising since as we showed in Paper I the stream velocity 
deprive the low mass halos of gas at high redraft. Therefore, this induce a large dependency on the behavior of the 
gas fraction in the low mass tail. To capture the low mass tail there is a range of intr oducing large erro rs in the gas 
fraction estimations, since these halos have less then 500 particles (as was shown in iNaoz et"aITl2009t ). Therefore, 
we introduce a wight function, th at account for the errors in estimating the gas fraction as a function of the number 
of particle in a halo, Nh- We use INaoz et al.l (120091 ) resolution study, their figure 4, and assume a scatter of 20% for 
halos that have more then 500 particles. For smaller nu mber of part i cles p er halo we adopted a simply linear function 
of the error of the gas fraction as can be estimated from INaoz et al.l (|2009( ). figure 4. Thus the wight function has the 
following form: 

wr/V^-/0-2 ifiV„> 500, 

- 1 7.5 X lO-^iV, + 0.575 if iV, < 500 . ^^^^ 

In figure [7| we show an example for the gas fraction as a function of mass. We depicted the gas fraction and the two 
fits as in Figure [H However, here we also show all of the points we considered in evaluating the fit (i.e., Nh > 100). 
The different values achieved for the two models are shown in Figure [H As depicted in this Figure, the N = 256 set 
produces a systematically low value for M,. even for the fbc = case. As mentioned in Paper I, this run suffer from 
low statistic, particularly in the large mass tail,. Furthermore, there was no convergence of the gas fraction, these 
are the main cause for this systematics. In addition as can be seen from Figure [51 Mc for the different values of fbc 
converge over the different ranges of redshift for the = 512 and A^ — 768 sets. 




Fig. 7. — A representative example of the gas fraction as a function of mass for the N = 512 set for three different redshifts (top to 
bottom) z = 15, 19 and z = 20. We compare between the new fit model, solid lines, using eq. ((2]|, and the old one, dotted line, using eq. i fTI 
We show all of the halo of which Nf^ > 100 (grey points) as well as the binned data points for A*^;, > 500. We consider the various values 
of the stream velocity (from left to right), v^,^ = 0, i;bc = lo'vbci ^bc = l-7(Tvbci a^r"! ''be = 3.4(Tvbc- 

The errors shown in Figures [2H4] and |8] are the maximum 1 — a errors between the two fit models. In addition we 
have used bootstrap method for the new model, in some cases of the = 512 set to test our evaluation of the error^, 
and found that they are consistent with choosing the maximum 1 — cr from the two fit models. We show the complete 
best fit parameters for the 512 set in table [21 they arc similar for the other sets. Note that for high rcdshift and large 
Wbc values the best fit parameters arc poorly constrain. Specifically the parameters f3 and 7 from equation ([2]) are 
sometimes so poorly constrains (i.e., more then an order of magnitude) that we omit the errors from the table, thus 
the symbol " " in table [2] means error larger then an order of magnitude. 

Note that in our calculation of the fit wc used /b,o which is the gas fraction in the high mass tail. This value is 
lower the the men cosmic baryonic fraction /b mainly because there might still be baryons which are bound_to the 
halo which arc missed sue to our halo finder algorithm (i )2.3|) . As can be seen in Figure [9] assuming /b.o ^ /b, grey 
lines, by either of the models results in a worse fit then using the calculated /b,o. 

N = 512 set for three different redshifts (top to bottom) z = 15, 19 and z — 20. We compare between the new fit 
model, using eq. ([2]), and the old one, using eq. ([TJ We show all of the halo of which Nh > 100 (grey points) as well 
as the binned data points for Nh > 500. We consider the various values of the stream velocity (from left to right). 
Vhc = 0, Vhc = iCTvbc, Vhc = l.TcTvbc, and Wbc = 3.4crvbc. 



We could not use the bootstrap method to all of the runs since .r 1, uii tii- 1 

^ m tew runs wc have somewhat low statistic sample. 
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Fig. 8. — Comparison between the new fit model and tiie old one. We consider in black the old model, using eq. {T}, and the color points 
are the fit for the new model, using eq. l(2j, calculated with minimum of 100 particles per halo (see text). Wc consider (from left to right) 
the N = 256,512 and N = 768 sets. We consider the various values of the stream velocity d^c = 0, Hbc = lo'vbci '''be = l-7o"vbci a-nd 
''be = 3.4(Tvbc (blue squares, purple diamonds, brown stars and red triangles respectively). 
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Fig. 9. — Comparison between the new fit model and the old one while using /b q and /]-,. We consider the N = 512 run for z = 15 
representative example. We compare between using value, solid (dotted) grey lines for the new (old) model. We also show the fit while 
using the /b,o value, solid (dotted) color lines for the new (old) model. We consider the various values of the stream velocity Hbc = 0, 
V\^c = lo'vbci ''be = l-7o"vbci s-nd Dbc = 3.4a'vbc (blue squares, purple diamonds, brown stars and red triangles respectively). 
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TABLE 2 

The best-fit parameters from equation ([JJ for the N = 512 set. 
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